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TECHNICAL NOTE 3572 

AMPLITUDE OF SUPERSONIC DIFFUSER FLOW PULSATIONS^ 
By William H. Sterbentz and Joseph Davids 


SUMMARY 

A theoretical method for evaluating the stability characteristics 
and the amplitude and the frequency of pulsation of ram- jet engines with- 
out heat addition is presented herein- Experimental verification of the 
theoretical results are included where data were available. 

Theory and experiment show that the pulsation amplitude of a high 
mass-flow-ratio diffuser Increases with decreasing mass flow. The theo- 
retical trends for changes in amplitude, frequency, and mean pressure 
recovery with changes in plenunt- chamber volume were experimentally 
confirmed . 

For perforated convergent-divergent- type diffusers, a stability 
hysteresis loop was predicted on the pressure- recovery, mass -flow-ratio 
curve. At a given mean mass-flow ratio, the higher value of mean pres- 
sure recovery corresponded to oscillatory flow in the diffuser while the 
lower branch was stable. This hysteresis has been observed experimentally. 

The theory indicates that for a ram- jet engine of given diameter, 
the an5)litude of pulsation of a supersonic diffuser is increased by 
decreasing the relative size of the plenum chamber with respect to the 
diffuser volume down to a critical value at which oscillations cease. In 
the region of these critical values, the stable mass-flow range of the 
diffuser may be increased either by decreasing the combustion- chamber 
volume or by increasing the length of the diffuser. 


INTRODUCTION 

Supersonic diffusers are usually selected so that at design mass 
flows pulsation- free operation will be obtained. However, there are con- 
ditions of operation, such as acceleration of maneuvering, which may 
require reduced mass flows in the unstable region. It is therefore impor- 
tant that both the conditions under which flow instability will be encoun- 
tered and the resulting pulsation amplitudes and frequency be known. 

^Supersedes NACA RM I52I24, "Amplitude of Supersonic Diffuser Flow 
Pulsations" by William H. Sterbentz and Joseph Davids, 1952. 
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Certain resonance criterions which reasonably predicted both the 
conditions for the occurrence of diffuser flow pulsations and the fre- 
quency of pulsation were evaluated assuming that a ram- jet engine acts, 
in effect, as a Helmholtz resonator (ref- l) . It was found that a dif- 
fuser would resonate if the subcritical diffuser characteristic had a 
positive slope of sufficient magnitude. At conditions where the slope 
of this cui*ve was less than the critical value, the ram- jet engine would 
not resonate. A similar criterion for the occurrence of diffuser pulsa- 
tion was suggested from analyses presented in reference 2. 

The occurrence of flow pulsation in external- compression diffusers 
is related in re^rence 3 to the velocity discontinuity or vortex sheet 
downstream of the intersection of two shock waves. Flow pulsations 
start Tdien this vortex sheet just enters or leaves the diffuser cowl- 
ing. The entrance of the vortex sheet into the diffuser cowling marks 
a change in the slope of the diffuser characteristic -vdoich arises from 
the entrance of increasing percentages of low-energy mass flows. The 
amount of change in diffuser characteristic required to induce flow 
pulsations in any ram— jet engine may be readily determined from the 
theoretical criterions of reference 1. 

Althou^ certain theoretical criterions as to the nature of dif- 
fuser pulsation have been evaluated, no theoretical evaluation of the 
variables affecting amplitude has been developed. This report, pre- 
pared at the MCA Lewis laboratory, presents a theoretical evaluation 
of the variables -viiich govern pulsation amplitude for the cold-flow 
engine. The theoretical method employed is a refinement of the 
diffuser-pulsation theory presented in reference 1. The method of 
reference 1 assumed small amplitudes of flow disturbances in order 
that quantities such as mass flow and the slope of the diffuser char- 
acteristic could be considered constant. 'When large finite amplitudes 
are allowed, these quantities must be recomputed in a point-to-point 
calculation throughout a cycle of oscillation. Several manerical com- 
putations which illustrate the nonlinear solutions for pulsation ampli- 
tudes and frequencies as a function of diffuser mass-flow ratio and 
geometry are included herein. Experimental data are presented for 
comparison with the theoretical trends where possible. 


SYMBOLS 

The following symbols and definitions were used in this report: 
A area, sq ft 

a velocity of sound, ft/sec 

f frequency, cycles/sec 
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M Mach number 

m mass flov, slugs/sec 

mo diffuser maximum mass flow, slugs/sec 

Am mass-flow perturbation, slugs/sec 

P total pressure, Ih/sq. ft abs 

AP pressure perturbation, Ib/sq ft 

APj pressure perturbation arising from the inertia properties of the 

perturbed flow, Ib/sq ft 

APjj pressure perturbation arising from changes in diffuser-shock 
position, Ib/sq ft 

t time, sec 

u velocity, ft/sec 

Au velocity perturbation, ft/sec 

V effective storage volume of ram jet " V^), cu ft 

effective air column, volume, cu ft 
total volume of ram jet, cu ft 
mass of virtual oscillating air column, slugs 
Wj, mass of gas occupying effective storage volume of ram jet, slugs 
y ratio of specific heats 

p density, slugs/cu ft 

Subscripts ; 
m mean 

0 eq[ulvalent free-stream condition 

1 diffuser inlet 

2 diffuser outlet 

4 exhaust -nozzle outlet 
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THEORETICAL AMLYSIS 

The "basic concept employed in the analysis of the amplitude of 
pulsation assumes that the resonant properties of a ram -jet engine are 
controlled "by the balance between the mass-flow-delivery characteris- 
tics of the diffuser, the storage capacity of the combustion or plenum 
chamber, and the discharge characteristics of the exhaust nozzle. For 
the engine under steady-flow conditions, the mass flow delivered to the 
combustion chamber is eq^ual to the mass flow discharged from the cham- 
ber. The presence of any transient disturbance induces an unbalance 
in the diffuser delivery and exhaust-nozzle discharge rates leading to 
a periodic storage and evacuation of air mass in the combustion chanber 
and concomitant cyclic pressixre oscillations . 

A linear differential equation of motion is developed, on the 
basis of the Helmholtz resonator concept (ref. l), through a system- 
atic calculation of the terms involved in the conservation-of-mass 
equation as functions of the velocity pertTarbation of the diffuser air 
column and mean-flow parameters of the engine. For conditions of 
unstable operation, this linear vibratory system predicts a continuous 
increase of pulsation to infinite amplitude. In actual systems, of 
course, a positive damping force becomes important at sufficiently 
large amplitudes, thus limiting the' oscillation. Such a positive damp- 
ing force arises from the nonlinear characteristics of the pressure 
changes incurred from shock oscillations during unsteady operation of 
the diffuser . 

From the conservation-of-mass equation 

Ang = Am^ + (l) 

the differential equation of motion developed in reference 1 is as 
follows : 

^c^A d^ ^ ^c- 1^2 r-1 /m\ „"1 d , 

dt2 " rPgagMslamg 2 [m^J A^ dt ^ + 

If second-order effects of terms invol'wlng Mg are neglected, the 
mass flow m^ may be stated as 
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“x = 


rPaA^Mg 

^2 


Miltlplying ty this expresslori for and nondimenslonallzlng ^^ 2/^2 

transforms equation (2) to 



r-1 

2y A^^Pg 


d 

dt ^ + 



Small oscillations were assumed in reference 1 so that the coefficients 
of the differential equation could he treated as constant values. 'When 
large finite amplitudes are allowed, the variation in quantities such as 

d(P2/Po) 


d^nig/i 


moJ 


} Pg, Sp, and m throughout the cycle must he included. How- 


ever, equation (2a) can still he applied in a point-to-point calculation 
assuming constant coefficients provided that the increments in the cycle 
are small enough. (The coefficients are, of course, the instantaneous 
values and thus are recomputed for each point.) 


To obtain results for specific diffusers requires a choice of a 
"frozen- state" pressure-recovery mass-flow-ratio curve. The "frozen- 
state" pressure recovery is defined as that recovery which would he 
obtained across the diffuser if steady-flow conditions could be imposed. 
This curve may generally be estimated from the specific diffuser geom- 
etry, the shock configuration, and the diffuser friction losses. 


For a hi^-mass-f low-ratio single-cone diffuser, peak frozen-state 
pressure recovery id.ll occur at critical mass flow. Pressure losses at 
this condition for such a diffuser will be Incurred from a shock con- 
figuration consisting of one oblique shock followed by a normal shock 
plus the estimated subsonic diffuser loss. For lower mass flows, the 
frozen-state pressure recovery decreases because decrements in mass flow 
from the critical condition are marked by the entrance of Increasing 
percentages of low-energy mass flows into the inlet as the vortex sheet 
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moves progressively closer io the diffueer centerhody. At zero mass 
flow, the diffuser pressure recovery will he approximately the recovery 
computed across a normal shock at free-stream Mach number. For the 
purposes of this report, it is considered sufficiently accurate to 
define the frozen-state diffuser pressxire recovery as varying linearly 
with mass -flow ratio from critical to zero mass flow for the high-mass - 
flow-ratio single-cone diffuser (fig. 1(a)). 

For low-mass -flow- ratio single-cone diffusers, peak pressure recov- 
ery will occur at some value of subcrltical mass -flow ratio. The range 
of mass -flow ratios from critical to peak pressure recovery is charac- 
terized by approximately constant total -pressure flow. At further 
reduced mass flows, the frozen-state pressiire recovery continuously 
decreases as the voirbex sheet enters the inlet and progressively advances 
toward the diffuser cone. By the method presented in reference 4, the 
entrance of the vortex sheet into the inlet and the decrement in pres- 
sure recovery with mass-flow ratio is easily confuted for any single- 
cone diffuser. The low-mass-flow curve can again be approximated by 
strai^t line segments as in figure l(b) . 

Certain cone-type diffusers encounter flow separation from the 
cone surface at reduced mass flows. At present it is not possible to 
predict "the mass-flow ratio at which this separation will occur or the 
decrement in pressure recovery it will induce. The method of refer- 
ence 5, however, does predict -vdiether or not a cone-type diffuser will 
encounter flow separation. 

Another type of frozen-state diffuser pressure-recovery mass-flow- 
ratio curve easily established is that curve commorly obtained with 
open-nosed, convergent-divergent, perforated diffusers. For this inlet, 
peak pressure recovery will occur at approximately critical mass flow. 

The theoretical frozen^state pressure recovery varies linearly along a 
line with decrements in mass-flow ratio from critical such as to main- 
tain constant Mach numbef at the diffuser throat until normal-shock 
pressure recovery at free-stream conditions is reached. This portion 
of the curve (fig. 1(c)) represents the shock in the convergent section 
of the diffuser. The pressure recovery then remains constant at normal- 
shock values with further decrements in mass-flow ratio and this portion 
of the curve represents the normal, shock ahead of the inlet. 

After the frozen-state diffuser characteristics are determined, the 
shape of the composite curve of motion and the mean diffuser character- 
istics under oscillating conditions are obtained in the following 
manner: In equation (2a) the following substitutions may be made: 

d 

Au = ^Au 


dAu . . dr 


dAu 


Aii = 


Au 


dt 
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Also^ the foUowing quantities are assumed to he equivalent: 


Wc = PgV 


% = P2^A 
2 ^^2 


Thus 


equation 


(2a) becomes , 


after simplification. 



Equation (3) describes a system of curves which relates the variation 
of air-column perturhation velocity Au with the air-column perturha- 
tion acceleration Au. This relation permits the drawing of tangential 
directions in the coordinate system for, the coordinates Au and Au. 
For example, for Au = 0, the slope dAu/dAu is infinite in the suh- 
critical flow region and for Au = 0, the slope dAu/dAu is equal to 
the quantity 



The proper values of pressure, tenperature, and mass flow occur- 
ring in the coefficients of equation (3) are determined hy successive 
apprxjxlmatlon at each point. The mass -flow change given hy equa- 
tion (7) of reference 1 stated in terms of Au and Au is 


Am 

“5 


1 

^^2 



(APi + APji) 




TPgAl 

™0®l 


Au 


(4) 


■where the pressure changes are 
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AP 


I 


4 


All 


or stating terms of P2^m^ APjj 



( 5 ) 


and APj-]- is eq^ual to the change in frozen- state pressure recovery 
from the mean. The change in a^ arising from adiahatic temperature 
changes is 



( 6 ) 


where &2 m and Pg axe the respective mean values of the diffuser- 
outlet acoustic velocity and pressure. 

An exaiiple of the use of equations (4), (5), and (6) follows. It 
is desired to establish the critical mass-flow curve on a ,Ail against 
An chart for some fixed mean value of mass-flow ratio. The origin on 
this chart, of course, represents the mean mass-flow conditions of the 
problem. The quantities immediately known are Am/mQ, Pg j^, 
ag jj^, the geometrical constants, and the variable Au. Quantities to 
be determined are APj, Pg, ag, and the appropriate values of the var- 
iable Au. 

By successive approximations using equations (5) and (6), the val- 
ues of APj, Pg, and ag are easily established for any fixed value 
of Au. (it should be noted that pressure Pg is the sum of the mean 
pressure Pg and the pressure changes Z^j and ZiPjj.) These quan- 
tities, when substituted in equation (4), permit the evaluation of Au 
and the establishment of one point on the Aii against Au chart rep- 
resenting critical flow conditions. When this calculation is repeated 
for other values of Au, the conqilete critical mass -flow curve is estab- 
lished. Other subciltical mass -flow-ratio curves are established in 
this same manner. 


r^ 

to 

in 

CV] 



2537 


NACA TN 3572 


9 


The eq^uation defining the supercritical mass-flow curve is 



P^moA^Au - 



= 0 


( 7 ) 


This equation is obtained from a substitution of equations (7), (ll)^ 
and (12) into equation (lO) of reference 1. For any given value of Au, 
the value of the Instantaneous pressure Pg and the acoustic veloc- 
ity ag may readily be determined by successive approximations from 
equation (4) . Thus, if the values of Pg and ag for a given value 
of Au are known, the proper value of Au may be computed from equa- 
tion ( 7 ) . The repetition of this calculation for other values of Au 
establishes the complete supercritical mass-flow-ratio curve. 

An illustrative solution of equation (3) is plotted in figure 2 
which presents a coordinate systan of Aa as a function of Au. The 
graph is divided into the subcritical and supercritical flow regions 
by the method just described. Starting at Au = - 600 feet per second 
and Au =0, the curve goes around in a clockwise direction until super- 
critical flow conditions in the diffuser are established,- and there- 
after describes a closed figure continuously. Also, starting from 
rest (Au = 0 and Zii = O), the same closed curve is reached after a 
number of increasing clockwise spiral turns. 

The variation of Au as a function of time t is obtained from 
a second graphical integration and the results of the exaiiple are shown 
in figure 3. The abscissa of a point on figure 2 is represented on the 
ordinate of figure 3 and the ordinate of a point on figure 2 is the 
slope of the curve on figure 3, Starting from the origin, the ampli- 
tude of Au gradually Increases with time until the value of the final 
or natural amplitude is reached. For the case where the initial dis- 
turbance is larger than the final value of amplitude, the pulsation 
diminishes with time to the final anplitude. By use of equations (4), 
( 5 ), and (6) and the graphical variation of Au and Au with time, 
the variation of ram- jet plenum-chamber pressure may be readily 
obtained. 

Some of the effects of changes in engine geometry, flow conditions, 
and diffuser characteristics win be evaluated and compared with avail- 
able experimental data in the following discussion. All theoretical 
amplitude and frequency curves presented were evaluated by the foregoing 
method. 
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EESDLTS AHD DISCUSSION 
Hi^-Mass -Flow-Ratio Diffuser 

CompEirisons of tlie theoretical and experimental pulsation ampli- 
tude and frequency for a Mgh-mass-flow- ratio diffuser with two plenum- 
chamber volumes at a free-stream Mach number of 2.0 is presented in 
figure 4. Schematic drawings of the ram- Jet engine and diffuser are 
shown in figure 5. One experimental frequency data point was used to 
establish the proper value of required in equation (3) and this 

point is tailed in figure 4(b) . The remainder of the frequency curves 
and the resulting pulsation amplitudes for both plenum- chamber volumes 
were then conputed using this value for the volume V^. The value of 
the ratio of the virtual oscillating air column to the total ram- Jet 
engine volume Tot the small plenum- chamber engine was 0.740 and 

for the large plenum^- chamber engine, 0.323. 

The theoretical variation and msignitude of pulsation anplltude 
with "both volumes is in good agreement with experiment over the range 
of mass flows investigated. The anplitude increased steadily with 
decreases in mass -flow ratio. In hoth cases, experimental frequency 
decreased gradually with increasing mass-flow ratio compared with the 
nearly constant value predicted hy theory. The experimental amplitude 
and frequency decreased as the plenum-chamber volume was increased as 
predicted by the theoiy. 

From the integration of the pressure-time diagram of a flow pulsa- 
tion (such as shown on fig. 6), the average increment or decrement in 
the frozen-state value of diffuser pressure recovery "which results from 
the pulsation can be conputed. The resulting curves are Included in 
figure 4(c) . Flow pulsation caused a decrease in pressure recovery 
below that value assumed for the frozen-state pressure recovery to a 
mass-flow ratio of about 0.61 for the small volume engine and "to about 
0.48 for the large volume engine, after "which an Increase was computed. 

A theoretical effect of an increase in diffuser pressure recovery caused 
by & decrease in plenum-chamber volume was also noted experimentally. 
However, "the experimental rise in pressure recovery "with mass-flow ratio 
in the range 0.46 < mg/niQ < 0.62 was not predicted by the theory. The 

discrepancies may be partly due to the crudity of the assumed frozen- 
state pressure-recovery mass-flow curve. 

No flow pulsations "were found by computation to be sustained for 
supercritical flow conditions eind all transient disturbances were found 
"to be danped out. 

The theoretical pressure-time "wave form only approximately predicts 
the experimental pressiare-time fluctuations as can be noted from fig- 
ure 6. Some variations from the experimental "values are to be expected 
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■because of the simplifying assumptions included in the analysis . A 
closer agreement of the theoretical and experimental wave forms at 
eq.uilihrlum hy a different analysis is presented in reference 6. 


Open-Nosed, Convergent-Divergent, Perforated Diffuser 

In the foregoing analysis of the properties of a high-mass-flow- 
ratio diffuser, small point-to-point increments in mass flow, pressure, 
and temperature were used to recompute the point-to-point values of the 
coefficients of equation (3) in obtaining the theoretical solutions. 

To consenre time and labor in the following solutions, only a few large 
increments of mass flow, pressure, and temperature were used. Thus, the 
general trends will still be the same although, if smaller increments 
were employred in the calculation, absolute values of the results would 
be somewhat changed. 

A theoretical solution of the pulsation characteristics of an 
arbitrary, hypothetical, open-nosed, convergent-divergent, perforated 
diffuser was obtained (see fig. 7). Two modes of operation were found 
to be theoretically possible. In the first mode, if the diffuser were 
started in the zero mass -flow condition with no flow oscillations, 
pulsation-free flow with approximately free-stream normal- shock recov- 
ery would persist over the mass -flow-ratio range to 0.50. At a mass- 
flow ratio of 0.50, however, flow pulsations will set in and the dif- 
fuser pressure recovery will increase several percentage points. Flow 
pulsations, once started, were found theoretically to persist over the 
entire subcrltical mass-flow range. This stability nysteresi's is 
defined, by theory, to arise from a critical pulsation amplitude such 
that flow perturbations less than this critical value are damped oxrt, 
•vdiereas perturbations greater than this critical value ^e amplified 
as an oscillatory disturbance to a new natural amplitude. For this 
particular hypothetical diffuser, the critical amplitude is approxi- 
mately the flow disturbance necessary to induce instantaneous values 
of the mass -flow ratio greater than 0.50. The variation of values of 
the critical and the natural amplitude with mass -flow ratio are shown 
on figure 7. 

The existence of a stability hysteresis and the rise in diffuser 
pressure recovery ^d-th the onset of flow pulsation are confirmed by 
experiments reported in reference 7. Reference 7 presents the perform- 
ance of a series of perforated, convergent-divergent diffusers having 
numerous contraction ratios and perforation distributions. Typical 
of the results of this investigation are the data plotted in figure 8 
for the 1.63 -e-contraction-ratio diffuser obtained from figure 17 (u) 
of reference 7. The lower branch of the hysteresis curve was generally 
obtained idien the tunnel flow was first Initiated. 'When the mass-flow 
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ratio was increased into tlie unstable region (from about 0.48 to 0.82)^ 
tbe mean-pressure recovery jumped to its oscillatory value. The dif- 
fuser then continued to operate on the upper or oscillatory curve even 
to nearly zero mass flow throu^ this diffuser. 


Some Additional Theoretical Trends 

The effects on -amplitude and frequency of varying plenum- chamber 
volume V or diffuser-air-column volume are sho"wn in -figures 9 

arifl 10. In figure 9, the plenum-chamber volume is -varied while the 
diffuser volume is held constant. As the plenum-chamber volume is pro- 
gressively decreased, the pulsation amplitude and frequency gradually 
Increase until a cidtlcal value is reached at -vdiich pulsations cease. 

In this region, the stable mass-flow range of a supersonic diffuser may 
be improved by further reductions in plenum- chamber volume. 

In figure 10, the euir-column or diffioser volume is varied -while 
the plenum-chamber volume is held constant. In this case, as the dif- 
fuser volume is Increased, the pulsation amplitude Increases and fre- 
quency decreases . Again, a critical ratio of plenum-chamber volume to 
diffuser volume exists such that further Increases in diffuser volume 
progressively increase the stable mass -flow-ratio range of the diffuser . 

• Thus, from figures 9 and 10, it can be noted that for any ram- jet 
engine of given diameter, the amplitude of pulsation of a supersonic 
diffuser may be decreased by increasing the relative size of the plenum- 
chamber volume -vd-th respect to the air-column or diffuser volume. Also, 
the frequency of pulsation for any ram-jet engine of given diameter will 
be decreased as either diffuser volume or plenum- chamber volume or 
both are increased. 


SUMMARY OF RESULTS 

The following summary presents the theoretical and experimental 
results of the effects of engine geometry and diffuser-flow conditions 
on supersonic diffuser pulsation amplitude and frequency. 

1. Theory and experiment show that pulsation amplitude of a high- 
mass-flow ratio diffuser increases -with decreasing mass-flow ratio. 

The theoretical trends for changes in amplitude, frequency, and^mean 
pressure recovery -vd.th changes in combustion- chamber -volume were 
experimentally confirmed. 

2. For perforated convergent-divergent-type diffusers a stability 
hysteresis loop was predicted on the pressure-recqvery miass -flow-ratio 
curve. At a given mean mass-flow ratio, the hl^6r value of mean-pressure 


2537 



2537 


MCA TIT 3572 


13 


recovery corresponded to oscillatory flow in the diffuser while the 
lower branch was stable. This hysteresis has been observed experimentally. 

3. The theory indicates that for a ram- jet engine of given diam- 
eter, the amplitude of pulsation of a supersonic diffuser is increased 
by decreasing the relative size of the plenum chamber with respect to 
the diffuser volume to a critical value at which oscillations cease. 

In the region of these critical values, the stable mass-flow range of 
the diffuser may be increased either by decreasing the combustion- 
chamber volume or by increasing the length of the diffuser. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 20, 1952 
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0 .2 .4 .6 .8 1.0 

Mass-flow ratio, 1112 /niQ 


(c) Open-nosed, convergent- divergent, perforated diffuser. 

Figure 1. - Typical variation of frozen-state pressure recovery 
with mass-flow ratio for three diffusers. 
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Figure 2. - Illustrative exan 5 >le of theoretical -growth of minute transient 

disturbance and diminution of large transient disturbance to natural amplitude 
of arbitrary ram Jet at fixed mean mass-flow ratio and free-stream Mach number. 
First integration. 



Berturbatlon velocity, Au, ft /sec 
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0 .1 .2 .3 .4 .5 .6 

Time, t, sec 

Figure 3. - Second integration of growth and diminution of trEinsient 
disturbances presented in figure 2. 
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Time, t, mllUsec 

(t) Mass-flow ratio: experiment, 0.41; theory, 0.40. 

Figure 6. - Comparison of theoretical pressure-time variation with 
experiment for hi^-mass-f low-ratio diffuser (small volume engine) . 
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Figure 8. - Pulsation hysteresis and effect of flow 
pulsation on diffuser pressure recovery. Perforated, 
convergent- divergent diffuser. Experimental inlet 1.65-e 
(ref . 7 ) . 
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(a) Frequency. 



(b) Amplitude. 

Figure 9. - Theoretical effect of plenum-chamber volume on diffuser stability, 
pulsation amplitude and frequency. Arbitrary, hypothetical ram-Jet engine 
with constant alr-colvnnn volume. 
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(a) Frequency. 







